Two distinct ancient components in the Sculptor Dwarf 
Spheroidal Galaxy: First Results from DART 1 
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ABSTRACT 

We have found evidence for the presence of two distinct ancient stellar com- 
ponents (both > 10 Gyr old) in the Sculptor dwarf spheroidal galaxy. We used 
the ESO Wide Field Imager (WFI) in conjunction with the VLT/FLAMES spec- 
trograph to study the properties of the resolved stellar population of Sculptor 



out to and beyond the tidal radius. We find that two components are discernible 
in the spatial distribution of Horizontal Branch stars in our imaging, and in the 
[Fe/H] and Vhei distributions for our large sample of spectroscopic measurements. 
They can be generally described as a "metal-poor" component ([Fe/H]< —1.7) 
and a "metal-rich" component ([Fe/H]> —1.7). The metal-poor stars are more 
spatially extended than the metal-rich stars, and they also appear to be kine- 
matically distinct. These results provide an important insight into the formation 
processes of small systems in the early universe and the conditions found there. 



Even this simplest of galaxies appears to have had a surprisingly complex early 



evolution. 
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Subject headings: dwarf galaxies: general — dwarf galaxies: individual(Sculptor 
Dwarf Spheroidal) 

1. Introduction 

The Sculptor (Scl) Dwarf Spheroidal galaxy (dSph) is a close companion of the Milky- 
Way at high galactic latitude (6 = —83°). Located at a distance of 72 ±5 kpc (Kunkel & 
Demers 1977), with v^ei — 109.9 ±1.4 km/s (Queloz, Dubath & Pasquini 1995) it has a low 
total (dynamical) mass, 1.4 ±0.6 x 1O 7 M (Queloz et al. 1995). Irwin & Hatzidimitriou 
(1995) made the first modern comprehensive determination of the physical properties of Scl. 
They determined a core radius, r c = 5.8 ±1.6 arcmin, and a tidal radius, r t = 76 ± 5 arcmin 
with a position angle of 99 ± 1 degrees, and confirmed it to be a moderately luminous 
galaxy in the dSph class (My = —10.7 ± 0.5), with a modest central surface brightness 
(X ,v = 23.5 ±0.5 mag/arcsec 2 ). Like most dSph galaxies Scl does not appear to have much, 
if any, HI gas (e.g., Bouchard, Carignan & Mashchenko 2003). 

Observations of Scl have revealed a sizeable population of RR Lyrae variable stars (e.g., 
Kaluzny et al. 1995), clearly indicating that its stellar population contains a globular cluster 
age component. Scl also shows tantalizing evidence in its extended horizontal branch (HB) 
for an unusual enrichment history (e.g., Hurley-Keller, Mateo & Grebel 1999; Majewski et 
al. 1999). Monkiewicz et al. (1999) have detected stars several magnitudes below the oldest 
possible main-sequence turn-offs. Their accurate photometry allowed them to conclude that 
the mean age of Scl is similar to that of a globular cluster, but that there is an age spread 
within this epoch indicative of an extended star formation history. Thus the entire star 
formation history of Scl apparently lasted only a few Gyr, and after this initial period of 
activity more than 10 Gyr ago Scl appears to have been dormant. 

There have been several previous spectroscopic studies of individual stars in Scl, firstly 
to determine the kinematic properties (e.g., Armandroff & Da Costa 1986; Aaronson & 
Olszewski 1987; Queloz et al. 1995) which have indicated a moderate mass-to-light ratio 
M/L= 9 ± 6(M/Ly) . This value is significantly higher than that typical for a globular 
cluster, but much lower than found for many dSphs (e.g., Draco and Ursa Minor). There 
have also been studies of the abundances of small samples of individual Red Giant Branch 
(RGB) stars at low and high resolution (Shetrone et al. 2003; Tolstoy et al. 2001, 2003). These 
stars were selected from the central region of Scl, and although they showed a wider range in 
[Fe/H] than would be naively expected based on colour-magnitude diagram (CMD) analysis, 
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no conclusions could be drawn on the large scale spatial variations of these properties. 

Here we present the first results from the DART (Dwarf Abundances and Radial veloci- 
ties Team) large programme at ESO. We present the results for \ hd and [Fe/H] measurements 
from our FLAMES spectroscopy of 401 RGB stars for the first galaxy in our sample, Scl 
dSph. The relatively high signal/noise, S/N (ps 10-20 per pixel) of our data has enabled us 
to derive both accurate metallicites (~ 0.1 dex from internal errors) and radial velocities (to 
±2 km/s). This is the first time that a large sample of both types of measurement has 
been made in a single galaxy. A more detailed description of our data and data reduction 
techniques will be provided in subsequent papers (Irwin et al. in prep; Hill et al. in prep; 
Battaglia et al. in prep). Here we provide an overview of our first intriguing results. 

2. Observations 
2.1. Imaging 

The ESO Wide Field Imager (WFI) observations were obtained from the ESO archive 
for the central region of Scl, and more extended observations of outer fields were made by us. 
We used the pipeline reduction software developed by the Cambridge Astronomical Survey 
Unit for processing mosaic camera imaging (Irwin & Lewis 2001; Irwin et al. 2004). These 
observations were used to study the photometric properties of Scl stars out to the tidal radius 
and as target selection for VLT/FLAMES observations of stars along the RGB. 

2.2. Spectroscopy 

All VLT/FLAMES observations were made in Medusa mode, giving ~120 fibres to be 
placed over a 25 arcmin diameter field of view. Most of the observations were taken with 
the Low Resolution grating LR8 (which includes the Ca II triplet (CaT) lines), apart from 
the central pointing where High Resolution gratings were used (see Hill et al. in prep). The 
CaT method is well calibrated to provide an accurate estimate of the [Fe/H] abundance of 
an RGB star (Rutledge, Hesser & Stetson 1997; Cole et al. 2004) and also to provide an 
accurate radial velocity measurement. The FLAMES fields were placed at varying distances 
from the centre of the galaxy to beyond the tidal radius. 

The FLAMES data were all reduced, extracted and wavelength calibrated using the 
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GIRBLDRS pipeline 1 provided by the FLAMES consortium (Geneva Observatory, Blecha et 
al. 2003). For sky subtraction, velocity and equivalent width (EW) estimation, we developed 
our own software which was thoroughly checked on multiple observations of the same fields 
taken at different times. Wavelength shifts were generally found to be negligible (ps ±1 
km/s). Each target spectrum was automatically continuum-corrected and cross-correlated 
with a zero-continuum Gaussian model CaT template. With a velocity determined it is then 
straightforward to estimate the EW of the individual CaT lines. We found that for reliable 
[Fe/H] determinations from CaT calibrations a S/N > 10 is required, and for reliable v he i a 
S/N > 7.5 is sufficient. All the plots in this paper use the S/N > 10 data set, as we always 
consider both the [Fe/H] and the v he i measurements. 

We were able to verify the zero point accuracy for both velocity and abundance deter- 
minations by comparison with independent previous measurements of stars in Scl, which we 
had deliberately reobserved. We also have overlap between our own HR and LR samples to 
verify consistency of velocity and [Fe/H] measurements using different gratings and based 
on different abundance determination methods. 



3. Analysis 

3.1. Spatial Distribution 

Previous studies have already suggested that the spatial distribution of the Horizontal 
Branch stars in Scl shows signs of a population gradient, with the red horizontal branch 
stars (RHB) being more tightly concentrated than the blue horizontal branch (BHB) stars 
(e.g., Hurley-Keller et al. 1999; Majewski et al. 1999). Unlike previous wide field studies, 
our WFI imaging data extends beyond the nominal tidal radius and with an average 5-cr 
limiting magnitude of V=23.5 and 1=22.5 also probes well below the horizontal branch. This 
has enabled us to unequivocally demonstrate that the BHB and RHB stars have markedly 
different spatial distributions (see Figure 1). 

These RHB and BHB populations were selected from the WFI CMD shown in Figure 2 
using the selection boxes marked. It is clear that the BHB region is almost completely 
uncontaminated by foreground stars, whereas the RHB zone includes a significant underlying 
foreground component. To make a quantitative estimate of the contamination of the RHB 
sample by foreground stars we selected a region of similar CMD foreground occupancy, well 
away from the Scl stellar locus (also outlined in Figure 2). The radial density distributions 
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of the RHB component were then corrected by the expected foreground contamination and 
the results for both the RHB and BHB distributions together with their density ratio as a 
function of radial distance, are shown in the lower panels of Figure 1. The different spatial 
occupancy of the two populations is striking and provides strong evidence that we are seeing 
two distinct components. The characteristics of the BHB and RHB are also consistent with 
different ages (e.g., <2 Gyr), or different metallicities (A[Fe/H]~0.7 dex), from theoretical 
modelling of globular cluster horizontal branches by Lee et al. (2001). 

3.2. Metallicity Distribution 

Our VLT/FLAMES spectroscopic sample provides measurements of [Fe/H] for 401 in- 
dividual candidate RGB-selected stars with final S/N>10, of which 308 have a high mem- 
bership probability (these are plotted in Figure 3), defined to be 80 km/s < Vhei < 150 km/s 
i.e., approximately within 3a of v sys (see Figure 4). At the high Galactic latitude of Scl \hei 
is a good estimator of membership probability assuming v sys = 110 km/s, and a ~ 10 km/s. 

For those RGB stars which were determined to have a high probability of membership, 
we show the distribution of [Fe/H] measurements as a function of elliptical radius (the 
equivalent distance along the semi-major axis) from the centre of Scl in the upper panel 
of Figure 3. A well-defined metallicity gradient is apparent with a similar scale size to the 
RHB versus BHB spatial distributions. The lower panel of this figure highlights the abrupt 
change in the [Fe/H] distribution that occurs at about r= 0.2 deg (2r c ) and also the broad 
range in [Fe/H] over the whole galaxy. Although both inner and outer distributions cover a 
broad range of values, the inner is centred at [Fe/H]= —1.4 with a range of ±0.6dex and the 
outer distribution is centred at [Fe/H] = —2 with a similar range of ±0.6dex. The drop in 
<[Fe/H]> at 2r c from the centre also corresponds to the transition between the dominance 
of the RHB and BHB populations seen in Figure 1. This is not seen if the colour of the RGB 
alone is taken as an indication of metallicity. 

3.3. Kinematics 

In the upper panel of Figure 4 we plot the sample of Vhei measurements versus elliptical 
radius for which the continuum S/N > 10. The distinction between foreground stars and 
members of Scl is generally unambiguous. The clustering around the systemic velocity of 
v sys = 110 km/s is clear and is highlighted by the dashed line and the dotted 3a deviation 
limits. Although present, foreground contamination can be seen to be almost negligible, 
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even in the outer regions, which simplifies subsequent analysis. 

We have divided our "member" sample into metal-poor ([Fe/H< —1.7) and metal-rich 
([Fe/H]> —1.7) distributions. This difference is quantified in the lower panels of Figure 4 
which shows the Vhei distributions over three radial zones for the different metallicity com- 
ponents separately. The variation in v sys and a in these three regions is computed using the 
maximum likelihood approach described in Hargreaves et al. (1994), but modified to directly 
compute the likelihood function over a 2D search grid. We found that v sys does not vary 
significantly between different populations or regions, but a does (see Figure 4). 

One obvious feature of the velocity distributions is the difference in the kinematic prop- 
erties of the two metallicity components. In general the more metal poor component has 
a higher velocity dispersion than the more metal rich component. This difference is most 
striking in the central panel (0.2° < r < 0.5°) of the triptych where a is 11. Oil. km/s 
and 6.8±0.8 km/s respectively. The swap in the dominance of the components is stark in 
the outer zone (r > 0.5), where there are hardly any stars at all in the metal-rich category 
and provides compelling evidence that we are indeed observing two dynamically different 
components in Scl. 

4. Conclusions 

Our FLAMES results show that Scl contains two distinct stellar components, one metal- 
rich, -0.9 > [Fe/H] > -1.7, and one metal-poor, -1.7 > [Fe/H] > -2.8. The metal-rich 
population is more centrally concentrated than the metal poor, and on average appears to 
have a lower velocity dispersion, a meta i-rich = 7± 1 km/s, whereas <J me tai- P oor = 11 ± 1 km/s 
(see Figure 4). Our WFI results show that there is a comparable difference over the same 
scale in the spatial distribution of BHB stars and RHB stars. Plausibly combining these 
results enables us to distinguish these two distinct components on the basis of kinematics 
(Figure 4), metallicity (Figure 3) and spatial distribution (Figure 1). These large samples 
allow us for the first time to quantify these large scale trends in a statistically meaningful 
fashion. 

There are indications that this is a common feature of dSph galaxies. Our prelimi- 
nary analysis of Vhei and [Fe/H] measurements in the other galaxies in our sample (Fornax 
and Sextans dSph, Battaglia et al., in prep) also shows very similar characteristics to Scl, 
especially in the most metal poor component. The mix of populations between the inner 
and outer regions is clearly changing. The radial velocity studies of Wilkinson et al. (2004) 
and Kleyna et al. (2004) have also considered the possibility that kinematically distinct 
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components exist in Ursa Minor, Draco and Sextans. 

What mechanism could create two ancient stellar components in a small dwarf spheroidal 
galaxy? A simple possibility is that the formation of these dSph galaxies began with an 
initial burst of star formation, resulting in a stellar population with a mean [Fe/H] < —2. 
Subsequent supernovae explosions from this initial episode could have been sufficient to cause 
gas (and metal) loss such that star formation was inhibited until the remaining gas could sink 
deeper into the centre and begin star formation again (e.g., Mori, Ferrara & Madau 2002). 
Thus the subsequent generation(s) of stars would inhabit a region closer to the centre of the 
galaxy, and have a higher average metallicity (> —1.4) and different kinematics. Another 
possible cause is external influences, such as minor mergers, or accretion of additional gas. 
It might also be that events surrounding the Epoch of Reionisation, strongly influenced the 
evolution of these small galaxies (e.g., Bullock, Kravtsov & Weinberg 2001) and resulted 
in the stripping or photoevaporation of the outer layers of gas in the dSph, meaning that 
subsequent more metal enhanced star formation occured only in the central regions. 

The full abundance analysis of the FLAMES HR data (Hill et al. in prep) will provide 
more details of the chemical enrichment history of Scl. This will hopefully enable us to 
distinguish between two episodes of star formation, each of which created the separate kine- 
matic components seen here, or more continuous star formation, manifested as a gradient in 
velocity dispersion and metallicity from the centre of the galaxy. 

Here we have shown that even dSphs, the smallest galactic systems, did not form simply, 
as might be assumed. The wealth of data presented here points to the presence of dynamically 
distinct components in these galactic building blocks, or templates, and provides a stimulus 
to theories of the formation and evolution of dSphs. 
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of Arts and Sciences, and thanks Renzo Sancisi for a careful reading of the draft. KAV thanks 
the NSF for support through a CAREER award, AST-9984073. MDS would like to thank 
the NSF for partial support under AST-0306884. 
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Fig. 1. — The distribution of Horizontal Branch stars from WFI imaging of the Scl dSph. The 
upper panels show the different spatial distributions of blue horizontal branch stars (BHB), 
and red horizontal branch stars (RHB) as selected from the M„, V— I Colour-Magnitude 
Diagram (CMD) shown in Figure 2. Also shown, to illustrate the foreground contamination 
in the RHB distribution, are a CMD-selected sample of foreground stars to match the RHB 
contamination density (see Figure 2). The ellipse is the tidal radius of Scl as defined by 
Irwin & Hatzidimitriou (1995). The lower panels quantitatively show the spatial variation of 
number density per unit area of these three population components (BHB in blue, RHB in 
red and foreground in green) and the ratio of BHB to RHB number densities after correcting 
the RHB component for foreground contamination. The physical scale at the distance of 
Sculptor is 1 degree to 1.4 kpc. 



This preprint was prepared with the AAS IATgX macros v5.2. 
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Fig. 2. — The Colour-Magnitude Diagram for the WFI coverage of Scl shown in Figure 1. 
The green and magenta symbols are the Red Giant Branch stars which were observed with 
VLT/FLAMES and for which we have accurate Vhei and [Fe/H] measurements (S/N>10). 
The potential members of Scl are shown as magenta dots and non-members shown as green 
stars. Also shown are the regions used to define the BHB, RHB and foreground comparison 
sample (FG). 
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Fig. 3. — VLT/FLAMES spectroscopic measurements of [Fe/H] for the 308 potential RGB 
velocity members of Scl versus radius (top panel). The lower panel shows the distribution of 
[Fe/H] for: the entire sample of Scl member RGB stars (dotted line); the 97 stars within the 
central, r=0.2 degree region (red solid line); and the 211 stars beyond r>0.2 degrees (blue 
solid line). 



-12- 



200 - 



-100 



VLT/FLAMES Scl dSph 



0.5 1 

Elliptical Radius (degrees) 



1.5 



< r < 0.2 



0.2 < r < 0.5 



0.5 < r < 1.8 



o=7.6±0.6 
0=9.3*1.2 



A 



ct=6.8±0.8 
o-=11 ±0.7 



mi 



= 14.6±4.0 
=8.S±1.0 



1 00 1 20 1 40 

Vhel 



1 00 1 20 1 40 

Vhel 



100 120 140 

Vhel 



Fig. 4. — The top panel shows Vhel as a function of elliptical radius for all stars satisfying S/N 
> 10. The 308 stars which are potential members are plotted as red stars ([Fe/H] > -1.7) 
and blue circles ([Fe/H] < -1.7), while the 93 green crosses are assumed to be non-members. 
The lower panels quantify how the velocity distribution changes as a function of metallicity 
and radius. The low metallicity distribution is plotted in blue, and the higher metallicity in 
red, as in the upper panel. 



